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Abstract 

Steady-state and transient computer models of the RL10A-3-3A rocket engine have been 
created using the ROCket Engine Transient Simulation (ROCETS) code. These models 
were created for several purposes. The RL10 engine is a critical component of past, 
present, and future space missions; the model will give NASA an in-house capability to 
simulate the performance of the engine under various operating conditions and mission 
profiles. The RL10 simulation activity is also an opportunity to further validate the 
ROCETS program. The ROCETS code is an important tool for modeling rocket engine 
systems at NASA Lewis. ROCETS provides, a modular and general framework for 
simulating the steady-state and transient behavior of any desired propulsion system. 
Although the ROCETS code is being used in a number of different analysis and design 
projects within NASA, it has not been extensively validated for any system using actual 
test data. The RL10A-3-3A has a ten year history of test and flight applications; it 
should provide sufficient data to validate the ROCETS program capability. The 
ROCETS models of the RL10 system were created using design information provided by 
Pratt &Whitney, the engine manufacturer. These models are in the process of being 
validated using test-stand and flight data. This paper includes a brief description of the 
models and comparison of preliminary simulation output against flight and test-stand 
data. 


Introduction 

The ability to simulate rocket engine system behavior is 
critical to the successful design and operation of such 
systems. Rocket engine cycles typically involve the 
interaction of several non-linear processes including 
combustion, heat transfer, flow resistance, turbine and 
pump operation. Relatively subtle variations of such 
processes can have a dramatic effect on the engine efficiency 
and stability. Predicting the behavior of rocket engines 
requires the use of digital computers and a great deal of 
effort on the part of scientists and programmers. 1 

The RL10A-3-3A rocket engine, manufactured by Pratt & 
Whitney, is an important component of commercial and 
military launches performed by the United States. The 
RL10 is used to propel the Centaur upper stage vehicles 
which carry payloads from a high altitude into orbit. The 


ability to adequately predict the behavior of the RL 10 has a 
significant impact on the success of such missions. The 
NASA Lewis Research Center has begun an effort to develop 
a model of the RL10 for in-house use. This model has been 
created using design data such as performance maps, line 
volumes and resistances, and heat transfer correlations 
supplied by Pratt & Whitney under contract to NASA (see 
also Ref.2). Ultimately, the goals of this effort include the 
development of more detailed RL10 component models using 
analysis tools which are still under development at NASA. 

NASA Lewis has created its version of the RL10A-3-3A 
computer model using the Rocket Engine Transient 
Simulator (ROCETS) software. 3 ROCETS has been selected 
for this task because it provides a modular, flexible, and 
powerful modeling capability in a non-proprietary package. 
The ROCETS program has also been used to model several 
other rocket engine systems, such as the Space Shuttle Main 
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Engine (SSME) and Space Transportation Main Engine 
(STME). The RL10 modeling activity is an opportunity 
to evaluate ROCETS capability to simulate the operation 
of expander cycle engines. Although the RL10 is not 
highly instrumented, it has a rich database of test-stand and 
flight operation which can be used to validate ROCETS 
model output 

ROCETS was designed to perform both steady-state and 
transient analysis. The first, steady-state version of the 
RL10 model includes volume, line, and rotor dynamics, but 
uses a steady-state heat transfer subroutine based on test 
data. This model was used to predict the performance of the 
engine at the full power level, and for several different fuel- 
to-oxidizer mixture ratios. These predictions have been 
compared with Pratt & Whitney simulations and test-stand 
data. 

The transient version of the RL10 model is derived from 
the steady-state version, but uses a five node transient heat 
transfer model suitable for engine-start simulations. The 
original combustion property tables have also been replaced 
with tables which cover the lower pressure region of 
operation (for start-up). The differences between the steady- 
state and transient model performance predictions at full 
power are discussed further in subsequent sections of this 
paper. The transient RL10 model has also been used to 
simulate the start behavior of the engine. Results from 
this simulation are compared with flight telemetry data. 

In the following sections of this paper, the RL10A-3-3A 
engine system and its model in ROCETS are discussed in 
greater detail. A comparison of steady-state model 
predictions with test-stand data is then presented. Finally, 
a comparison of start transient behavior with flight data is 
discussed 

Description of RL10A-3-3A Rocket Engine 

The RL10A-3-3 A is a cryogenic hydrogen/oxygen rocket 
engine manufactured by Pratt & Whitney for upper stage 
and orbital transfer applications. The engine system, 
represented schematically in Figure 1, is based on the 
expander cycle. In this cycle, heat from the main 
combustion chamber is transferred to the fuel via a 
regenerative cooling jacket; the heated fluid is used to drive 
a turbine. Power from the turbine drives the fuel and 
oxidizer pumps, connected through a gearbox. The engine 
develops a rated thrust of 16500 lbf (73395 N), with a 
specific impulse of approximately 445 seconds. The fuel 
turbopump has a design speed of 30000 rpm while the 
LOX pump design speed is 12000 rpm. Combustion 
chamber pressure is nominally 475 psia. The engine 
consumes approximately 6 Ibm/sec of hydrogen and 30 
Ibm/sec of oxygen at rated thrust. 2 

In its present configuration, the RL10A-3-3 A is not 


intended to be throttled. Chamber pressure is controlled by 
varying the amount of turbine bypass flow; this is achieved 
through a pneumatic servo-system which prevents the system 
from overshooting its desired operating point on start-up. 
Mixture ratio control is possible in order to optimize 
propellant utilization during the mission. During start, the 
engine accelerates to full power using only tank pressure and 
heat picked up in the cooling jacket; no boost pumps are 
currently used. 

Description of RL10 Model in ROCETS 

The ROCETS code is a general-purpose system analysis 
code. This code is written in ANSI standard FORTRAN 77 
and provides a flexible, modular, and powerful system 
modeling capability. The ROCETS package includes a 
library of commonly used component models such as 
dynamic volumes, flow resistances, pumps, turbines, etc. as 
well as physical property tables for hydrogen, oxygen, 
helium, nitrogen, and combustion products. The user may 
also attach new generic or special purpose component 
subroutines. The model configuration is defined in a text 
format, which is then automatically converted into a 
FORTRAN program. A separate set of instructions directs 
execution of the model, allowing the same model 
configuration to be used in a number of different 
applications, including parametric design, steady-state 
analysis, transient simulation, and generation of linearized 
models. 3 The ROCETS system dramatically reduces the time 
required to generate and debug model simulation programs. 

The ROCETS model of the RL10A-3-3A system is depicted 
in Figure 2. Each major component is represented by a 
separate subroutine in this model. Duct volumes and 
resistances, however, are lumped such that each resistance or 
volume in the model may actually represent multiple 
sections of ducts or manifolds. The volumes and resistances 
are lumped according to recommendations made by Pratt & 
Whitney based on their own RL10 engine models. 

The first and second stage pump performances are represented 
by maps (tables) of head coefficient and efficiency versus 
flow coefficient 4 These maps are based on test-stand data. 

At present the maps do not contain data for the high and low 
extremes of flow-coefficient The LOX pump is represented 
as fifth-order polynomial curve-fit relationships of head 
coefficient and efficiency to flow coefficient Approximate 
leakage rates from the pump stages to the gearbox have been 
included. The fuel turbine performance is represented as a 
map of efficiency versus velocity ratio. 4 The turbine 
pressure drop is modeled in a fashion similar to compressible 
flow through an orifice. 

The steady-state and transient versions of the RL10 model 
use different heat transfer calculations. The steady-state 
version uses a map (table) of heat transfer rate versus 
chamber pressure and mixture ratio based on empirical test 


2 



data. The transient version of the model uses a more 
complex five-node heat transfer model. In the transient 
model, the hot-gas side heat transfer coefficients are based 
on a Bartz correlation 5 while the cool-side coefficients are 
determined using a Colburn correlation 6 . The thermal 
capacity of the metal is taken into account, as well as the 
volume dynamics of the cooling jacket. 

Valves and ducts are both modeled as fluid resistances. 
Effective areas can be varied to simulate the action of the 
valves. Actuator dynamics (forces required to actuate 
valves, pneumatic and hydraulic line pressures, etc.) have 
not been included in the current model. The flow area of 
the fuel turbine bypass valve (FTB V) is used to set 
chamber pressure and thrust. The mixture ratio control 
valve (MRV) is used to set mixture ratio. For steady-state 
analysis, it is possible to set these valve areas to a desired 
value or automatically vary the valve areas until a desired 
target thrust and mixture ratio are achieved. For start 
transient operation, the areas of all valves are set according 
to schedules rather than using closed loop control. 

Each dynamic volume module has two state variables: 
density and internal energy. The combustion chamber 
pressure, temperature, and mixture ratio, and the fuel pump 
shaft speed and discharge flow-rate are also system states. 
These states represent the key dynamic parameters in the 
system model. Several additional algebraic balances, or 
‘loops*, are also required for this system. The states and 
additional balances define the Jacobian matrix for the 
system. This matrix is inverted and used in an iterative 
Newton-Raphson algorithm to determine model solution. 
For steady-state analyses, the state time-derivatives are 
driven to zero. During transient simulations, the 
integration corrector-term error is driven to zero. 3 

Steady-state model predictions vs. Test data 

In order to verify the steady-state RL10A-3-3A model, 
system performance was predicted for several different 
mixture ratios about the rated thrust level of 16500 lbf. 
These steady-state results were compared with the Pratt & 
Whitney test-stand data. Tables 1 through 5 show the 
ROCETS model predictions for ten key performance 
parameters at five different thrust and mixture ratio settings. 
The degree of difference between the model and test data is 
also listed in the tables. Each parameter has been classified 
according to the percent difference found between test and 
model results (in 1% intervals). The total number of 
parameters in each percent-error group is then totaled over 
the five test cases. These percent totals are shown as a 
histogram in Figure 3. 

As indicated by Figure 3, the model predictions are within 
4 % of the test data for all measured performance parameters. 
The figure also indicates that the majority of these predicted 
values are within 1% of test results, averaged over the five 


test cases. 

The transient version of the RL10 model was also driven to 
steady-state for comparison with test data. The principle 
difference between the steady-state and transient models is in 
the cooling jacket heat transfer calculation. The steady-slate 
model heat transfer is based on test data while the transient 
model uses a more complex theoretical calculation as 
described in the introduction. The results from the transient 
model do not match the steady-state test data quite as well as 
the steady-state model does. The RL10 transient model still 
matches measured test data within 4%, with the exception of 
turbine inlet temperature, which matches within 7% in all 
tests. 

Start Transient simulation results y^Flighl 
Data 

To add transient simulation capabilities to the model used for 
steady-state analysis above, it was necessary only to replace 
the steady-state heat transfer model with a dynamic heat 
transfer model and to extend the hot-gas property tables to 
include lower chamber pressures; the original version of the 
model already included volume dynamics, and fluid and rotor 
inertias. Because the basic RL10A-3-3A engine is not 
normally throttled, there is no suitable empirical data for 
mainstage transient operation. In order to validate the 
ROCETS RL10 model for transient operation, it is 
necessary, therefore, to compare simulation results with data 
for start-up and shut-down operations. Thus far, efforts have 
been focused on modeling the start sequence only. 

Simulating the start transient operation of an engine is 
extremely difficult for three reasons: 1) the system pressures 
are low enough for two-phase flow to exist, 2) the low speed 
performances of the RL 10-3-3 A turbines and pumps are not 
well characterized, and 3) the near zero values of system 
states leads to numerical instabilities in the simulation itself. 

Figures 4, 5,6, and 7 show plots of chamber pressure, 
venturi inlet pressure, LOX pump speed, and LOX pump 
discharge pressure respectively. On each plot, simulation 
results are overlaid with data from Atlas/Centaur flight AC- 
69. In each of these plots, the ROCETS prediction parallels 
the flight data but the acceleration begins too late relative to 
flight measurements. The simulation results shown here are 
representative, although it has been observed that the start- 
time may be earlier or later depending on the choice of initial 
conditions and other factors as discussed below. 

Discussion of Results 

The close correlation of the RL10 steady-state model output 
with test data, as shown in Figure 3, indicates that ROCETS 
can accurately predict the performance of the system at 
mainstage operation. There is good agreement between the 
various system pressures, temperatures, flows and shaft 
speeds. It should be noted here that the turbine bypass 
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(FTB V) and mixture ratio (MRV) valves are varied, in the 
engine and in the model, in order to achieve the desired 
thrust and inlet mixture ratio for the engine. The areas 
themselves are not actually measured and so the correlation 
with test data is not known precisely. 

It has also been noted that using the transient heat transfer 
model to predict mainstage steady-state performance is not 
as accurate as using the heat transfer correlations derived 
from empirical data. Although the 7% disagreement in 
temperatures is not an uncommon degree of error for such 
models, efforts are currendy underway at NASA Lewis to 
improve the transient heat transfer model. 

RL10 modeling efforts at NASA Lewis are presently 
focused on refining the start transient model. The time at 
which the system begins accelerating to full-power 
operation is extremely sensitive to variations in the flow 
rate through the Interstage Cooldown Valve, and in the fuel 
pump discharge pressure. The three areas of the RL10 
model presendy being given closer scrutiny are 1) the 
turbopump performance maps, 2) the resistance at the 
pump discharge cutwaters (where the pump impeller and 
discharge volute meet), and 3) the flow model for the 
Cooldown Valves. For the first 0.30 seconds of time after 
the start signal, there is an appreciable flow through the 
pumps (0.5 to 1.5 Ibm/sec) while the shaft speed is still 
relatively low (less than 2000 rpm for the fuel pump). 

This condition causes the pump to actually dissipate energy 
from the liquid, creating a pressure drop across the pump. 
This area of operation has not been extensively tested nor 
documented for the RL10A-3-3A pumps, and the empirical 
pump maps do not describe this dissipative effect There 
may also be an additional fluid resistance due to cavitation 
in the pump discharge cutwaters. The cavitation losses are 
included in the model as simple resistances when the pump 
discharge pressures approach the fluid vapor pressures. The 
Cooldown Valves, particularly the Interstage valve (see 
Figure 1), play a major role in the start-transient The 
hydrogen will vaporize as it is vented overboard but the 
location where the vaporization occurs will vary depending 
on upstream conditions. These effects significantly 
complicate transient simulation of the start process. If the 
model predicts too high a pressure drop over the pump, 
there will be insufficient energy to bootstrap; if too little 
pressure drop is predicted, the turbine will receive too much 
power at first and the system will accelerate prematurely. 
The RL10A-3-3A engine is minimally instrumented, 
giving limited information with which to understand the 
true operation of the pumps at start-up. Numerical 
instabilities during the simulation have also thwarted 
efforts to determine the correct starting conditions and 
design parameters via parametric studies. Present RL10 
modeling efforts are focused on solving these problems. 

In addition to the acceleration delay, there are differences in 
the shapes of the simulated and flight data curves as the 


engine approaches full thrust. The RL10A-3-3A uses a 
pneumatic servo-control system to open the turbine bypass 
valve (FTB V) at higher chamber pressures. A model of this 
servo-control has not been available for inclusion in the 
model. It is believed that actuator dynamics are the cause of 
the differences in data as the system levels off to full power. 
The small mismatches in the steady-state levels at the end of 
the start transient occur because the start simulation is run 
without closed loop control of thrust level. 

Concluding R?mqrR s 

Steady-state and transient computer models of the RL10A-3- 
3 A have been created using the ROCETS software. This 
effort serves the double purpose of providing NASA with an 
in-house capability to predict the performance of the RL10 
engine, and of providing a test case to validate the 
capabilities of the ROCETS code. The models were created 
using design information provided to NASA by Pratt & 
Whitney. These models have been used to perform steady- 
state performance predictions and to simulate the system 
transient behavior during start-up. The predictions from the 
steady-state model are within 4 % of test data for all measured 
performance parameters, and most parameters are within 1% 
of the test data. The success of the steady-state analyses 
demonstrates the utility of this model for mainstage 
performance predictions. The start transient simulation 
results accurately depict the acceleration of the engine, but 
the time at which the acceleration process begins differs from 
that observed in flight Current efforts are focused on 
refining the component and system models, particularly in 
the area of the fuel turbopump, in order to improve the start 
predictions. With some additional refinements, the transient 
model will also prove useful for predicting the transient 
behavior of the engine. The ROCETS code has proven very 
useful in these modeling efforts, facilitating the creation of 
the RL10 model with minimal time and effort required to 
program and debug the simulations. 
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Tab!# 1 - T#*1 Point 1 - ROCETS vs T#st Dala 



ROCETS 

Test Point 1 

ROCETS vs Test 

Medal Input 
THRUST (Ibf) 

16603.00 



Inlet O/F - Mixture Ratio 

5.63 




27.00 




36.60 




35.80 




174.70 



Output 

Fuel Pump Disch Pressure (static psia) 


1061.40 

■■■MM 

Fuel Pump Inlet Volumetric Flow (gpm) 

586.87 

580.60 


Venturi Inlet Pressure (static psia) 

799,95 

800.45 


Turbine Inlet Temperature (deg R) 

425.09 

439.32 

*3.24% 

Fuel Injector Manifold Pressure (static psia) 

527.96 

534.19 


LOX Pump Discharge Pressure (static psia) 

598.18 

613.80 

•2.54% 

LOX Pump Shaft Speed (rpm) 

12578.00 

12697.00 

-0.94% 

LOX Pump Inlet Volumetric Flow (gpm) 

206.99 

205.82 

0.57% 

LOX Injector Plennum Pressure (static psia) 

525.90 

532.50 

-1.24% 

Combustion Pressure (Inj. face static psia) 

468.72 

472.35 

-0.77% 


Tab!# 2 - T#st Point 2 - ROCETS vs T#st Data 



ROCETS 

Test Point 2 

ROCETS vs Test 

Model Input 
THRUST (Ibf) • 

16588.00 



Inlet O/F * Mixture Ratio 

5.55 



FUEL TANK Pressure (static psia) 

27.00 



FUEL TANK Temperature (deg R) 

38.60 



LOX TANK Pressure (static psia) 

35.30 



LOX TANK Temperature (deg R) 

174.70 



Output 

Fuel Pump Disch Pressure (static psia) 

1039.00 

1068.90 

-2.80% 

Fuel Pump Inlet Volumetric Row (gpm) 

593.44 

588.90 

0.77% 

Venturi Inlet Pressure (static psia) 

803.95 

808.56 

-0.57% 

Turbine Inlet Temperature (deg R) 

419.37 

432.92 

-3.13% 

Fuel Injector Manifold Pressure (static psia) 

528.89 

535.67 

•1.27% 

LOX Pump Discharge Pressure (static psia) 

601.86 

619.10 

*2.78% 

LOX Pump Shaft Speed (rpm) 

12617.60 

12747.00 

-1.02% 

LOX Pump Inlet Volumetric Flow (gpm) 

206.08 

205.59 

0.24% 

LOX Injector Plennum Pressure (static psia) 

525.91 

532.67 

-1.27% 

Combustion Pressure (Inj. face static psia) 

469.23 

472.96 

-0.79% 


Tab!# 3 - T«st Point 3 - ROCETS vs Test Data 



ROCETS 

T#st Point 3 

ROCETS vs T#st 

Model Input 
THRUST (Ibf) 

16458.00 



Inlet O/F - Mixture Ratio 

4.99 



FUEL TANK Pressure (static psia) 

26.50 



FUEL TANK Temperature (deg R) 

38.50 



LOX TANK Pressure (static psia) 

35.80 



LOX TANK Temperature (deg R) 

174.70 



Output 

Fuel Pump Disch Pressure (static psia) 

1083.44 

1112.50 

1 

Fuel Pump Inlet Volumetric Row (gpm) 

640.56 

634.90 

: s 

Venturi Inlet Pressure (static psia) 

836.14 

839.53 


Turbine Net Temperature (deg R) 

382.29 

390.35 

•2.07% 

Fuel Injector Manifold Pressure (static psia) 

536.47 

540.20 


LOX Pump Discharge Pressure (static psia) 

634.20 

647.70 

•2.08% 

LOX Pump Shaft Speed (rpm) 

12921.20 

13040.10 


LOX Pump Inlet Volumetric Flow (gpm) 

200.47 

199.61 

0.43% 

LOX Injector Plennum Pressure (static psia) 

527.44 

531.50 

-0.76% 

Combustion Pressure (Inj. face static psia) 

473.75 

475.48 

-0.36% 
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Table 4 - Test Point 4 - ROCETS vs Test Data 


Modal Input 
THRUST (Ibf) 

Inlet O/F - Mixture Ratio 
FUEL TANK Pressure (static psia) 
FUEL TANK Temperature (deg R) 
LOX TANK Pressure (static psia) 
LOX TANK Temperature (deg R) 



Table 5 - Test Point 5 - ROCETS vs Test Data 



Model Input 
THRUST (Ibf) 

Inlet O/F - Mixture Ratio 
FUEL TANK Pressure (static psia) 

FUEL TANK Temperature (deg R) 

LOX TANK Pressure (static psia) 

LOX TANK Temperature (deg R) 

Output 

Fuel Pump Disch Pressure (static psia) 

Fuel Pump Inlet Volumetric Flow (gpm) 
Venturi Inlet Pressure (static psia) 

Turbine Inlet Temperature (deg R) 

Fuel Injector Manifold Pressure (static psia) 
LOX Pump Discharge Pressure (static psia) 
LOX Pump Shaft Speed (rpm) 

LOX Pump Inlet Volumetric Flow (gpm) 

LOX Injector Plennum Pressure (static psia) 
Combustion Pressure (lnj. face static psia) 


ROCETS 


Test Point 4 


ROCETS vs Test 
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Figure 1-RL10A-3-3A Simplified Schematic 
















Figure 2 - RL10A-3-3A model configuration in ROCETS 
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Figure 3 - Difference between ROCETS 
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tion provided by Pratt & Whitney, the engine manufacturer. These models are in the process of being validated using test- 
stand and flight data. This paper includes a brief description of the models and comparison of preliminary simulation 
output against flight and test-stand data. 
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